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ABSTRACT 

We study inhomogeneous magnetised cosmologies through the post-recombination 
era in the framework of Newtonian gravity and the ideal-magnetohydrodynamic limit. 
The nonlinear kinematic and dynamic equations are derived and linearised around the 
Newtonian counterpart of the Einstein-de Sitter universe. This allows for a direct com¬ 
parison with the earlier relativistic treatments of the issue. Focusing on the evolution 
of linear density perturbations, we provide new analytic solutions which include the 
effects of the magnetic pressure as well as those of the held’s tension. We confirm 
that the pressure of field inhibits the growth of density distortions and can induce a 
purely magnetic Jeans length. On scales larger than the aforementioned characteris¬ 
tic length the inhomogeneities grow, though slower than in non-magnetised universes. 
Wavelengths smaller than the magnetic Jeans length typically oscillate with decreasing 
amplitude. We also identify a narrow range of scales, just below the Jeans length, where 
the perturbations exhibit a slower power-law decay. In all cases, the effect of the field is 
proportional to its strength and increases as we move to progressively smaller lengths. 


1. Introduction 


The origin and the potential implications of the large-scale magnetic fields seen in the uni¬ 


particular, bypass these helds despite their widespread presence 

(Kronbere 

1994 

Garilli V Tavlor 

2002: 

Han V Wielebinskv 

2002: 

Beck 

2008: Valee J.P. 

2011) 

. The fact that the galactic, the 


cluster and the protogalactic magnetic helds have comparable strengths, between 10 ^ and 10 ® 
Gauss, could be inter preted as a sign that cosmic magnetism has cosmological origin, but the is¬ 


sue is far from c losed (jGrasso &: Rubinstein Il200ll : IWidrow II2002I : iKandus. Kunze &: Tsagas 112011 


Rvu et al Il2012h . Gosmological magnetic {B) fields could have played a role during both the early 
and the late stages of structure formation, since they introduce new ingredients to their environ¬ 
ment. During the dust epoch, for example, large-scale R-helds (if present) are essentially the sole 
source of pressure support. The latter is known to inhibit the growth of density perturbations and 
also determine the scale of the Jeans length, which in turn decides the size of the first structures to 
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collapse gravitationally. An additional, rather unique, magnetic property is their tension, reflecting 
the elasticity of the field lines and their tendency to react against any agent that distorts them. 


In the present article we analyse the effects of large-scale magnetic fields on the evolution of 
baryonic density inhomogeneities after recombination and until the moment the recent accelerated 
expansion of the universe starts. We do so by confining to subhorizon scales and use the Newtonian 
version of the 1+3 covariant approach to cosmological perturbation theory (jEllis &: van Elst 1119981 : 


Tsagas. Challinor &: Maartens I l2008l b After a brief introduction to the formalism, we proceed 


to analyse the magnetic component within the framework of the ideal magnetohydrodynamics 
(MHD). The latter is believed to provide an accurate approximation of the observable universe, 
which appe ars to be a very good e l ectrical conductor. We follow the steps of earlier relativistic 


treat ments ( Tsagas fc Barrow 1997 . 1998 : Tsagas fc Maartens Ibnnnl: ! Isarrow. Maartens fc Tsagas 


20071 ). which also allows for a direct comparison between the two approaches. At the same time. 


the relative simplicity of the Newtonian analysis enables a closer and more detailed look into 
the magnetic effects on density perturbations after recombination. We begin by providing the 
nonlinear expressions monitoring the evolution of inhomogeneities in the density distribution of 
the (highly conductive) baryonic matter. These are subsequently linearised around an Einstein-de 
Sitter background that also contains a weak and fully random magnetic field. Focusing on the role 
and the implications of the magnetic pressure, we confirm that it establishes a purely magnetic 
Jeans length, the scale of which depends on the field’s strength. We also provide a new (to the 
best of our knowledge) set of analytic solutions, revealing that the magnetic effect on linear density 
inhomogeneities is scale-dependent. More specifically, perturbations larger than the aforementioned 
Jeans length grow, but at a pace slower than in non-magnetised models. Wavelengths smaller than 
the Jeans length, on the other hand, typically oscillate with a time-decaying amplitude. Finally, 
we find that there is a narrow margin of scales, just inside the Jeans length, where the density 
inhomogeneities experience a (relatively) slow power-law decay. In all cases, the magnetic presence 
inhibits the growth of density perturbations, but the exact nature of the effect depends on the 
field’s strength and on the scale under consideration. 


Adopting the Newtonian framework also allowed us to incorporate the effects of the magnetic 
tension into the linear equations (see Appendix A). As we mentioned above, these effects carry 
the reaction of the field lines to distortions caused by changes in their shape and their orientation. 
Nevertheless, at the linear perturbative level, the role of the magnetic tension seems secondary com¬ 
pared to that of the field’s pressure. Finally, for completeness and also for comparison, we consider 
and discuss the magnetic effects on linear density perturbations assuming a pure (i.e. magnetic-free) 
Einstein-de Sitter background (see Appendix B). 


2. Newtonian covariant magnetohydrodynamics 


The covariant formalism 


see 


Ellis & van Elst 

(1998b 

Tsaeas. Challinor & Maartens 

(2998! 


for reviews) has been applied to cosmology by many authors. Electromagnetic fields have been stud- 
























- 3 - 


i ed us ing covariant techniques in lEhlers (|l96lh ; lEllis (|l973l i ; iTsaeas (j2005l i ; ISovrou Sz Tsaeas 


(200S) and the f ormalism was applied to magne t ised c osmological perturbations in lTsaeas &: Barrow 


19971 1 (see also lBarrow. Maartens fc Tsagas (|2007l l for a review). Nevertheless, a Newtonian co¬ 


variant study of perturbed magnetised cosmologies has been missing from the literature. 


2.1. Kinematics 


We use general coordinates {x“}, with a = 1,2,3, to dehne the metric tensor {hab) of the 
Euclidean space. The metric and its inverse matrix satisfy the constraints hach'^b = ^ab^ with 
6ab being the Kronecker symbol. In a Cartesian frame, where hab = the covariant and con- 
travaria nt components coi ncide and spatial derivatives are replaced by ordinary partial derivatives 
(e.g. see lEllis (jl97ll . Il973l ll. We adopt a fluid approximation, where the matter and the comoving 
(fundamental) observers are moving with velocity Va tangent to the flow lines. The ‘convective’ 
derivative of a given tensorial quantity T is therefore given by T = dtT + v'^daT, where the first 
term on the right-hand side describes changes at a fixed point and the second takes into account 
the motion of the fluid. Hence, the convective derivative of the fluid velocity, 

Va = dtVa + v’^dbVa , (1) 

provides the inertial acceleration. 


The complete kinematic description of the fluid, as well as that of the associated comoving 
observers, follows from the spatial gradient of the velocity field. Like any second rank tensor, the 
derivative dbVa decomposes as 

dbVa 2 ®^ab T '^ab T ^ab ) (2) 


where 0 = d^Va, (Jab = 5^6W> aiid ‘^ab = The scalar 0 describes variations in the fluid 

volume, while the tensors aab (shear) and ojab (vorticity) monitor changes in the shape (under 
constant volume) and in the orientation (i.e. rotation) respectively. Positive values of 0 indicate 
expansion and negative contraction. The volume scalar also dehnes a representative length scale 
(a) according to d/a = 0/3. In cosmological studies, the aforementioned length coincides with the 
scale factor of the universe. 


2.2. Hydrodynamics 


Within the Newtonian framework, the gravitational held is described by the potential (<h), 
which obeys the Poisson equation 


= - Kp — A. 


( 3 ) 


^Round brackets in the indices denote symmetrisation, square indicate antisymmetrisation and angled ones define 
the symmetric and trace-free part of second-rank tensors. Therefore, dihVa) = d(bVa) — {d'^Vc/A)5ab- 
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where is the Laplacian, p is the matter density and k = SttG. Also, the cosmological 

constant (A) is measured in units of inverse-time squared. Given that da^ is the gravitational 
acceleration and Va its inertial counterpart (see Eq. the vector 

Aa = Uq + 5(1$ , (4) 


carries the combined action of gravitational and i nertial forces. Th erefore, Aa is the Newtonian 


analogue of the relativistic 4-acceleration vector (jEllis Il97ll . Il97,‘ll ). The density and the total 


acceleration of the fluid satisfy their conservation laws, namely the continuity equation and the 
Navier-Stokes formula, which take the covariant form 

p = -Op and pAa = -daP - d’^iTab , (5) 


respectively. Here, p is the isotropic and Tiab is the anisotropic pressure of the medium (with 
TTafe = 7r{ab>)- Finally, one also needs an equation of state relating the pressure to the density. Since 
we will be dealing with non-relativistic matter, the pressure is thermal in nature and the associated 
sound speed is very small (compared to that of light, namely = dp/dp <C 1). 


To complete the hydrodynamic description, we need the propagation equations of the three 
kinematic variables and an equal number of c onstra i nts. F or o ur purposes, the key forrn ula is the 
propagation equation of the volume scalar (see lEllis (jl97ll ) and ISpvrou fc Tsagas (120081 ) for more 
details), namely the Newtonian version of Raychaudhuri’s equation. 


1 


1 


0 — —— 0 — — up -\- d'^Aa — 2((T — lo ) -|- A , 


( 6 ) 


with cT^ = cTafecr“*/2 and = 0Jab^°‘^/“^ giving the shear and the vorticity magnitudes respectively. 


2.3. Ideal magnetohydrodynamics 

In an electrically conducting medium, the charges are carried by the positive ions and by the 
electrons. Assuming global electrical neutrality an d adopting a single -f luid approach, the veloci ty 
of the matter is that of the centre of mass (e.g. see iGiovannini (120041 1: ISovrou h Tsagas (j2008l H 

1 


Va = 


(m+u+ -b m-Va ), 


( 7 ) 


m+ + rti- 

where m± are the masses of the ions and the electrons and their velocities. Then, at the ideal- 
MHD limit the electrical conductivity diverges and Ohm’s law reduces to a simple relation between 
the electric (Eq) and the magnetic field (Ea)- 


Ea = -eabcV^B^ , 


( 8 ) 


where eabc is the Levi- Civita tensor. On u sing the above, the magnetic induction equation takes 
the covariant form (see ISovrou Sz Tsagas (j2008l l for the details) 


Ba — —— QBa + {(Tab + (Vab)B 


( 9 ) 
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guaranteeing that the magnetic vector connects the same particles at all times (|Ellis Ill973l h In 
other words, when the ideal-MHD limit holds, the i?-field is frozen into the matter. 


The magn etic held is a source of isotropic and anisotropic pressure (e.g. see marker (jl979l l 


Mestel (Il999l )). with both contributing to the Lorentz force. A t the ideal-MHD limit, the latter 
recasts the Navier-Stokes equation into (ISpvrou fc Tsagas 1120081 ) 


pAa = -daP - d^TTab “ ^ daS"^ + B^’ShBa , (10) 

which replaces the standard expression given in Eq. ([ 5 ) 3 ). Note that the third term on the right- 
hand side is due to the magnetic pressure and the fourth comes from the held’s tension. When 
these two stresses balance each other out, the Lorentz force vanishes and the H-held can reach 
equilibrium. The effect of the tension term vanishes when the field lines are ‘geodesics’, in which 
case B^dbBa = 0. 


In the ideal-MHD limit, the matter and the H-held are conserved separately. In particular, 
the fluid continuity equation retains its ‘non-magnetised’ form (see Eq. ([5^) earlier). On the other 
hand, contracting the induction equation (see expression ([9])) along the held vector gives 


^B^y = -_eB^ -2aabll^^. 


( 11 ) 


where liah = (B^/ 3)6nh — BgBh- The abov e can be seen as the conservation law of the magnetic 
isotropic pressure (|SDvrou fc: Tsagas ll2008l L Similarly, starting again from Eq. Q, one obtains the 
conservation law of the held’s anisotropic pressure, namely 


4 2 o 

Ilaft = 0II(j6 -|- b) ~ 6 } ~ 2 H ^ab ■ 


( 12 ) 


3. Inhomogeneous magnetised cosmologies 

Magnetic helds appear to be everywhere in the universe. Galaxies and galactic clusters, in 
particular, support large-scale H-helds with strengths around 10“® G and 10“^ G respectively. 
Although the origin of cosmic magnetism is still a mystery, it is conceivable that H-helds could 
have played some role during the process of structure formation. 


3.1. The key variables 


The observed large-scale structure of the universe is believed to be the result of gravitational 
instability, a physical mechanism that allows small inhomogeneities in the density distribution 
of the matter to grow with time. In order to study the magnetic effects on the evolution of 
density perturbations, we need to deh ne a set of key var i ables that descri be these distortions 


Following the relativistic treatments of iTsagas &: Barrow (|l997l . Il998l l and iTsagas &: Maartens 
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([ 2 OOOI) , inhomogeneities in the density distribution of a magnetised medium are monitored through 
the dimensionless gradients 


Aa = - daP , 
P 


and 


Ba = —, daB^ 


(13) 


The former of the above describes spatial variations in the density of the baryonic matter, as 
measured between two neighbouring comoving observers (i.e. flow lines), while the latter does the 
same for the magnetic pressure. To close the system, one also needs the auxiliary variable 


Za = ada<d , 

describing spatial variations in the volume expansion. 


(14) 


3.2. The nonlinear eqnations 


The evolution of density inhomogeneities is determined by the convective derivative of A^. 
Starting from definition (fT3]l . using the continuity equation (see expression ([5^)), and then applying 
the commutation law between spatial gradients and convective derivatives, givej§ 


Aa — Z(j (o'ba T ^6a)A 


(16) 


Unlike the relativistic case (compare to Eq. (5.7.4) in iBarrow. Maartens &: Tsagas (|2007l B. matter 
perturbations are not coupled to those in the magnetic pressure at this stage. Instead, the coupling 
of the two occurs in the second-order differential equation for Aa- The reason lies in the treatment 
of space and time, which in the Newtonian theory are completely separate entities. 


Similarly, the nonlinear propagation formula of the expansion gradients comes from definition 
(d- The dot-derivative of the latter, together with Raychaudhuri’s equation (see Q) and the 
commutation law (fTSjl leads to 

Za = -^ 0Za - ^ KpAa + oAa “ {(Tba + ^^ba)2^ “ ‘2'ada{cr‘^ - U?) , (17) 

where Aa = da{d^Ai,) by dehnition. Although the above contains no explicit magnetic terms, the 
influence of the R-field is incorporated within the acceleration term Aa, as the Navier-Stokes 
equation ensures (see expression d). Note that the relativistic analysis allows for a direct 
coupling between the expansion and the magnetic pressure gradients (compare to Eq. (5.7.5) 


®By means of definition © and decomposition the commutation law between convective time-derivatives 
(i.e. “dot-derivatives”) and spatial gradients takes the form 

da^ = {datj))' + i Oda(j) + Aba + U]ba)d'’(j) . 

Note that the above (Newtonian) law also holds when ij) is replaced by an arbitrary vector or tensor field. 


( 15 ) 







- 7 - 


. There is no such coupling in (1171) (as well as in Eq. (|16p l 
and the reason is partly the aforementioned difference is the treatment of space and time between 
the two theories, and partly the absence of a magnetic contribution to the total gravitational mass 
of our Newtonian system. 

Finally, taking the convective derivative of (jl3b ). using the conservation law of the magnetic 
pressure (see Eq. (jlip l. the nonlinear commutation law (|15p and expression ()16p . we obtain 

Ba = + + (18) 


m 


Barrow. Maartens &: Tsaeas ( 20071 11 


The above monitors spatial variations in the (isotropic) magnetic pressure, between a pair of neigh¬ 
bouring observers. We point out the direct connection between the magnetic and the matter inho¬ 
mogeneities, reflected in the first two terms of Eq. (|18p . This coupling, which also occurs in relativis¬ 


tic stu dies - compare the above to the propagation formula (5.7.6) in iBarrow. Maartens fc Tsagas 
(120071 1. will prove particularly useful at the linear level. 


4. The linear regime 

The set (jl 6 p - (jl 8 p governs the nonlinear evolution of density inhomogeneities in a magnetised 
Newtonian medium of high electrical conductivity and zero total charge. In this section we will 
linearise these expressions around the Newtonian analogue of the Einstein-de Sitter universe. 


4.1. The unperturbed background 


Let us assume a homogeneous and isotropic background universe filled with a barotropic per¬ 
fect fluid (with p = p{p) and Tiab = 0 ) and allow for a weak and completely random ( i .e. sta tistically 
homogeneous and isotropic) magnetic field (see also iBarrow. Maartens fc Tsagas (120071 )). Later 
we will also consider a pure (i.e. a magnetic-free) Einstein-de Sitter unperturbed model (see Ap¬ 
pendix B). Assuming a random B-field means that Ba = 0 on average in the background. On the 
other hand, the magnetic pressure does not average to zero (i.e. 7 ^ 0 on average) and has a 

time dependence only. The randomness of the B-field guarantees that its presence does not destroy 
the uniformity of our zero-order model, although it adds to the total pressure of the system. In 
addition, the overall weakness of the field ensures that the Alfvfo speed (defined by = B"^ jp) 
is always well below unityEI Thus, ignoring the cosmological constant, our unperturbed model is 


^Demanding that cj 1 places no real constraint on the magnetic field. For instance, setting p ~ GeV^, 

which is the critical density of the universe today, translates into B <C 10““^ G at present. The latter condition is 
satisfied by all the known large-scale magnetic fields. 














governed by the se' 


p = -3Hp, H'^ = ^kp, H = -^kp (19) 

and 

[B’^y =, ( 20 ) 

since 0 = 3H = 3a/a to zero perturbative order. Solving ([Eb) for the scale factor, leads to 
a oc which is the familiar evolution law of the Einstein-de Sitter universe. Also, following 
(fEh i and (j20|) . gives p oc a~^ and B^ oc a~^, guaranteeing that the magnetic pressure decays faster 
than the density of the matter. Finally, it is straightforward to show that 

= (21) 

to zero order. The above implies that oc a~^ oc to ensure that the Alfven speed decreases 
inversely proportionally with the dimensions of the universe. 


4.2. Linear magnetised pertnrbations 


The background model determines our zero-order solution. Perturbatively speaking, quantities 
with nonzero background value will comprise our zero-order variables, while those vanishing there 
will be treated as first-order perturbations. When linearising, terms of perturbative order higher 
than the first will be neglected. Thus, given our background, the only zero-order quantities are the 
Hubble parameter {H), the density of the matter (p) and the magnetic pressure (H^). All the rest 
are of order one perturbatively. Note that differentiation (temporal or spatial) leaves the order of 
a perturbed quantity unaffected. 

Applying the above described linearisation scheme to our basic equations, namely to the non¬ 
linear formulae (fT6]l - (fT8]l . the latter reduce to 


and 


Aq — 2a 


Za = -2HZa - - UpAa aAa 


= 3 , 


( 22 ) 

(23) 

(24) 


respec tively. In contrast with the relativistic treatment (e.g. see § 6.3.4 in iBarrow. Maartens &: Tsaeas 


(|2007l i for comparison), the magnetic effects on linear density perturbations propagate only through 


^Hereafter, barred variables will indicate background quantities of zero perturbative order. 
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the acceleration term on the right-hand side of Eq. (1231) . Recalling that Aa = da{d^Ah), assuming 
an ideal fluid (i.e. setting iiab = 0 in (fT0]) i and using definitions (fT3]) we obtain 

aAa = - ]- cld^Ba +‘^da {cTs - , (25) 

2 p 

to linear order@ Note that = P/p ^ 1 dehnes the unperturbed adiabatic sound speed and 
/p (with <C 1 at all times) the background Alfen speed. Also, Ug = d(^i,Ba)d^^B°''>/2 and 
~ by construction. These variables may be seen as the magnetic analogues of 

the shear and of the vorticity respectively. Finally, following Eq. (1241) . the linearised gradients of 
the magnetic pressure evolve in tune with those in the density of the matter. More specifically 

l3a = ^Aa+C, (26) 

with C representing the integration constant. The system (I22n - (l24p monitors the linear evolution 
of density inhomogeneities in a weakly magnetised Newtonian medium of zero electrical resistivity. 
Alternatively, one can take the convective derivative of (l2^ and then employ Eqs. ([2^ . (f^ll to 
eliminate Za- Finally, using the auxiliary expression (I25h and the linear result (I26p . the aforemen¬ 
tioned set is replaced by 

Aa = -2HAa + ^ KpAa + + \ ~ ^ ~ 

which is a wave-like differential equation with additional terms due to the universal expansion, the 
presence of matter and the action of the magnetic field. 


4.3. The magnetic effects 


To linear order, the effects of the R-field on density perturbations propagate via the last two 
terms on the right-hand side of Eq. ()27ll . The former is due to the (positive) magnetic pressure 
and the latter is the result of the magnetic tension. Recall that the tension reflects the reaction 
of the field lines against any agent distorting their equilibrium. As a result, the tension effects 
are ex pected to become more pr ominent when the entanglement of the field lines is sufficiently 
large (jTashiro fc Sugivama 1120061 ). Put another way, although the magnetic shear and vorticity 
stresses (i.e. the terms cr^ = /2 and = d]f,Ba\d'^^B‘^^/2 respectively) are first-order 

perturbatively, they are expected to play a secondary role, unless the magnetic configuration is 
highly distorted. Thus, ignoring (momentarily only - see Appendix A) these two terms and then 


®Perturbatively, the gradient daB^ is treated as a first-order quantity, since it vanishes in the background. Then, 
given that daB^ — 2B^daBi,, both Ba and dtBa are of perturbative order half. This makes cr^ and lu% of order one. 
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harmonically decomposing the remaining expression give^ 


^(n) = + 7^/^P 


1-x (c^ + xc^ 1 ( ^ 


A 


(n) 


(28) 


where A// = 1/H is the Hubble radius and An = a/n is the scale of the inhomogeneity. Note 
that, we ha ve arrived to a homogeneous d ifferential equation, a s oppo sed to an inhomogeneous 
one (e.g. see Kim. Qlinto &: Rosner (1996); Tashiro fc Sugivama (2006))]^ This is also true when 
the tension stresses have been incorporated into the system (see Appendix A below). As a result, 
extracting analytical solutions becomes a relatively straightforward process, which is an attractive 
feature of the covariant formalism. 


It is worth pointing out that one can use Eq. (1281) even when the universe is dominated by 
pressureless non-baryonic cold dark matter (CDM), provided there is no relative velocity between 
the baryonic and the dark components. This happens because the magnetic field does not interact 
with the CDM, which contributes only to the total gravitational field through Raychaudhuri’s 
formula (see Eq. dH) in § 12.2p . Therefore, qualitatively speaking, the only change would be in 
the background density (i.e. in p - see expression (1281) 1. which would now correspond to the dark 
sector. Nevertheless, the background dynamics remain unaffected, since both the baryonic and the 
dark-matter densities obey the same evolution lawl§ On the other hand, the background dynamics 
will change drastically if we assume that the universe is dominated by dark energy or by a positive 
cosmological constant. Consequently, given that current observations support the idea of a recent 
universal acceleration, Eq. (1281) holds from the time of recombination (at z ~ 10^) until the start 
of the accelerated phase (at z ~ 0.5). 

An additional issue is the working domain of expression (I28p and the spectrum of the back¬ 
ground magnetic field. Since we use Newtonian gravity, we can only consider subhorizon-sized 
perturbations (i.e. those with A„ < Ah). Also, if one wants to focus on the magnetic effects and 
bypass those of the baryonic pressure, as we are going to do next, the scales of interest should not 
be very small and probably not much smaller than those of a galactic cluster. On all these lengths, 
the spectrum of the background (the random) magnetic field has been treated as scale-invariant. 

Expression (f28l) describes magnetosonic waves, propagating with an effective “sound” speed 
equal to + (2/3)c^ (see also Eq. (|T7)) ). Consequently, the R-field adds to the total pressure 


^We obtained (1281) by setting Aa = with representing standard vector harmonic functions 

and n > 0 the comoving wave numb e r of t h e associated mode. Also, by constru ction, daP^{n) = 0 = and 

(e.g. see Ellis ( 199flh : Tsagas. Challinor fc Maartens ( 2008l l). 


“A homogeneous differential equation, but of order four, was also obtained by iRuzmaikina fc Ruzmaikin (Il97ll l. 

®The role of CDM should become more prominent at the nonlinear level, when the baryonic matter starts falling 
into the ‘potential wells’ of the dark matter. T his increases the strength of the magnetic field, which remains frozen-i n 


with the baryons, beyond the adiabatic limit ( Bruni. Maartens fc Tsagas I I 2 OO 3 I : IPolag. Bartelmann fc Leach II 199^ . 
The enhanced _B-field could then backreact and have a stronger effect on the evolution of density perturbations. 
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of the system, which in turn increases the associated Jeans length. In fact, even when dealing with 
pressure-free matter, there is a purely magnetic Jeans length given by 


Aj = - Ca^H ■ 


(29) 


S 7.4.1 in 

Barrow. Maartens &: Tsaeas 

Kim. Olinto &: Rosner 

(1996lh 


Tashiro &: Sueivama ( 20061 ') 1. Note that Ca oc and \h oc which imply that Xj oc a oc 
In other words, the physical size of the magnetic Jeans length increases in tune with the dimen¬ 
sions of the universe. Consequently, when the first stars where formed, between z ~ 15 and 2 ; ~ 20 
(e.g. see ITashiro &: Sugivama (j2006l ll. the magnetic Jeans length (and the physical size of the first 
structures allowed to collapse gravitationally) was approximately 15 to 20 times smaller than it is 
today. Similarly, at the start of the accelerated phase (at 2 ; ~ 0.5) the magnetic Jeans length was 
slightly (roughly 1.5 times) smaller than its current value. 


Using definition (I29p one can estimate the scale of the magnetic Jeans length for some char¬ 
acteristic values of the background il-field. The same length-scale also provides the size of the 
first gravitationally-bound structures to form, in a scenario where the only pressure support comes 
from cosmological magnetic fields. For instance, setting B ~ 10“^ G, which is the typical strength 
of the cluster magnetic fields and also the maximum allowed by primordial nucleosynthesis, gives 
~ 10“®^ GeV^. Then, if we set the background density equal to the critical density (i.e. for 
P = Pc ^ 10“^'^ GeV"^, we find that Aj 1 Mpc. The latter is the size of a small galac¬ 
tic cluster. Alternatively, for B ~ 10“® G, which is the upper limit of a homogeneous B- 
field allowed by the Gosmic Microwave Background (GMB), the magnetic Jeans length drops to 
Aj 10 Kpc. Fi nally, recent surveys have suggested the presence of intergalactic magn etic fields 
close to 10“^® G ( Tavecchio et al 2010 : Ando fc Kusenko 201C : Neronov fc Vovk 201C 1. Adopt¬ 
ing this value, expression (|2^ gives Aj ~ 10“^ pc. Note that, as we have pointed out earlier, our 
analysis holds until the onset of the universal acceleration (at 2 ; ~ 0.5), at which moment the above 
given values are slightly smaller (by roughly 1.5 times - see earlier discussion). Overall, assuming 
that cosmological magnetic fields are the only sources of pressure support after recombination and 
given the presumed weakness of such field s, the first structures to form will have typical sizes of a 
small galactic cluster or smaller (see also iKim. Olinto &: Rosner (jl996l ll. In general, the B-field 
adds to the ambient thermal pressure of the matter and thus increases the effective sound-speed 
and the associated Jeans length by a small amount. 


Ignoring the matter pressure means that Eq. (I28p can be solved analytically. In particular, 
given that H = 2/3t and up = 4/3t^ (for both baryonic dust and CDM), the aforementioned 
differential equation reads 


dt^ 


4 dA(^) 
3t dt 



(30) 


The latter, which on subhorizon scales and in the absence of curvature effects agrees fully with its 
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relativistic analogue (e.g. see § 7.4.1 in iBarrow. Maartens &: Tsaeas (|2007l B. accepts the solution 


= Cl ts (725^-1) ^ ^-1 {C25^+l) ^ 


(31) 


with 



(32) 


being a positive constant (recall that = B^/p oc a~^ and Xn/^n oc a^/^). The above monitors 
the effects of cosmological B-fields on linear inhomogeneities in the density of baryonic dust. Note 
that, given the weakness of the magnetic fields (i.e. since = B‘^/p <C 1), the typical values of the 
a are smaller than unity. It is conceivable, however that on sufficiently small wavelengths (i.e. those 
with (A/i/An)^ > (75/32)c“^), the value of a may exceed unity (see third case next). Having said 
that, we should also stress that on very small scales the nonlinear effects start becoming important 
and our linear approximation is expected to break down. 


The first point to emphasise is that when there is no B-field (i.e. for a = 0), solution (1311) 
reduces to its non-magnetised (dust) counterpart, where +C 2 t~^ on all scales. Other¬ 

wise, the magnetic effect depends on the value of the a, namely on the strength of the field and on 
the scale of the perturbed mode (see definition (l3^ above). We may therefore distinguish between 
three alternative cases: 


• When 0 < a < 1, one can easily show that An > Aj and 0 < (\/25 — 24a — l)/6 < 2/3. 
Therefore, on scales larger than the magnetic Jeans length, solution (131 p reduces to a power- 
law that contains one growing mode. The latter, however, increases slower than in non- 
magnetised cosmologies. 


• Assuming that 1 < a < 25/24, namely that ^24/25Aj < An < Aj, gives —1/6 < (V25 — 24a— 
l)/6 < 0. Within this narrow margin of wavelengths, density perturbations either remain 
constant or decay at a pace no faster than A(n) oc It looks like a ’transition‘ range, 

where the gravitational pull is not strong enough to make the perturbations grow and the 
support of the field’s pressure is unable to make them oscillate (see next case). To the best 
of our knowledge, this type of behaviour has not been noted before. If observed, it could in 
principle provide evidence of a magnetic presence on relatively large scales. It is conceivable, 
however, that other sources of pressure could mimic this magnetic effect. 


• Finally, for a > 25/24 we have An < y^24/25Aj and 25 —24a < 0. Then, setting \/25^^'24a = 
^^/24a — 25 in the right-hand side of (j31h . shows that the inhomogeneity oscillates with de¬ 
creasing amplitude. More specifically, when a > 25/24, solution (13Ih takes the oscillatory 
form 


A(„) = 


Cl COS ^-\/24a — 25 In tJ -|-C 2 sin ^-\/24a — 25 In t 


(33) 


while its amplitude drops as As we mentioned above, solution (|33p holds on small scales 

(i.e. those with (A/i/A„)^ > (75/32)c“^ <(4> A„ < y^24/25Aj), where the magnetic pressure is 
strong enough to support against the gravitational pull of the matter. 
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In summary, one could argue that the -B-field has a negative effect on the growth of linear den¬ 
sity inhomogeneities, since the latter either grow slower or decay with time. Ana logous (negative) 
effects were observed both in relativistic studies fe .g. seelTsaeas &: Barrow (Il998li fo r the radiation 
era), as well as in Newtonian treatments (e.g. see iRuzmaikina &: Ruzmaikin (jl97l|)). In all cases, 
the reason is the extra pressure that the field introduces into the system and the overall magnetic 
effect is proportional to its strength. Although, the latter drops with time, the overall effect is 
decided by the a-parameter which is scale dependent (see definition ([32])) and remains constant. 


5. Discussion 

The continuous detection of large-scale magnetic fields in the universe, especially those found 
in galaxy clusters and in distant proto-galactic formations, adds support to the idea of primordial 
magnetism. Assuming that i?-fields have cosmological origin, it is conceivable that they could 
have played some role during the process of structure formation. In the present article we have 
re-considered the implications of large-scale magnetic fields for the linear evolution of density per¬ 
turbations, in the baryonic component of the matter, after recombination. This allowed us to 
employ a Newtonian approach, which is believed to provide an adequate treatment of the problem 
in hand on subhorizon scales and once the universe has gone past equipartition and decoupling. 
Technically, we proceeded by adopting the Newtonian version of the 1-1-3 covariant approach to 
cosmology, which also facilitated the direct comparison of our study with the previous relativistic 
treatments. In addition, one should in principle be able to exploit the techniques developed here, 
as well as the experience gained, to improve the general relativistic study of the issue. For instance, 
the technique used to incorporate the tension effects of the i?-field and still arrive at a homogeneous 
differential equation for the density perturbation (see Appendix A), is likely to work within the 
relativistic framework as well. 

We looked at the problem by adopting two alternative scenarios. The first allowed for a weak 
and fully random (statistically homogeneous and isotropic) H-field in our Einstein-de Sitter back¬ 
ground, while the second employed a pure (magnetic-free) unperturbed model. In the latter case 
the field did not have any practical effect on the linear evolution of density inhomogeneities, since 
it simply added a constant mode in between the growing and the decaying modes of the non- 
magnetised study (see Appendix B). Introducing a weak magnetic field into the background, on 
the other hand, strengthened its overall impact and lead to explicit magnetic terms in the solutions. 
The latter were found to differ qualitatively, depending on the strength of the B-field and on the 
scale of the perturbed mode. Roughly speaking, on wavelengths longer than the so-called magnetic 
Jeans length, density perturbations grow but at a pace slower than in the magnetic-free models. 
Inhomogeneities spanning scales smaller than the aforementioned Jeans length, however, oscillate 
with time-decreasing amplitude. There is also a narrow window of wavelengths, just inside the 
magnetic Jeans length, where the perturbations simply decay with time. Finally, incorporating 
the effects of the magnetic tension did not seem to cause any significant change to our results (see 
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Appendix A). Practically speaking, whether the aforementioned effects could leave an observable 
signature depends on the actual strength of the magnetic fields and on the sensitivity of the in¬ 
struments. Qualitatively speaking, we expect to see an overall ‘negative’ effect, since the field’s 
presence either slows down the standard growth-rate of the inhomogeneities, or forces them to 
decay. Analogous (negative) effects were also observed in the relativistic study of the radiation era 
and in all cases the reason is the extra pressure that the i?-held introduces into the system. 

It is also worth mentioning that our Newtonian analysis also confirmed the existence of a purely 
magnetic Jeans length, the scale of which depends on the strength of the held. This characteristic 
length-scale determines the size of the hrst structures to form gravitationally, assuming that mag¬ 
netic helds are the only sources of pressure support in the post-recombination universe. Given the 
presumed weakness of cosmological magnetic helds, we expect that the aforementioned hrst struc¬ 
tures are small, with their typical sizes never exceeding that of a small galactic cluster. Generally 
speaking, the U-held adds to the ambient thermal pressure of the matter and thus increases the 
total effective Jeans length by a small amount. 


A. Incorporating the magnetic tension 


The magnetic effects analysed in section § 14.31 account for the (positive) pressure of the held 
only. To incorporate the tension component of the Lorentz force (see decomposition (IlOp in S 12.3p . 
one needs to include the magnetic shear and the magnetic vorticity terms (a^ and respectively) 
on the right-hand side of Eq. (|27l) . Taking the time-derivative of the latter, using a set of aux¬ 
iliary relations (see Eqs. (IA3p and (IA4h below) and then harmonically decomposing the resulting 
expression, gives 


d^A 


(n) 


dt3 


10 d^A(„) 2(1 + a) dA(„) 

3t dt^ 


(Al) 


3t^ dt ’ 

where a = {XjjXnf' = constant (see dehnition The above solves to give 

A(,) = Cl d (v^55^-1) + (C25^+1) ^ ^ (A2 ) 

which is identical to solution ()3ip with the exception of an additional constant mode on its right- 
hand side. The latter carries the effects of the magnetic tension and plays a role only when 
a > 2b124, in which case the other two modes decay with time (see also § 14.31 above). Over¬ 
all, when d ealing with linear density per t urbat ions, the magnetic tension is typically a secondary 
player (see iBarrow. Maartens &: Tsaeas (j2007l l for a discussion of the relativistic case). On the 
other hand, the tension plays the prominent role when studying the magnetic effects on linear ro ¬ 
tational distortions, bo t h in t he Newtonian and the relativistic framework fsee IWasserman (jl978l k 
Dosopoulou fc Tsagas ( 2014 1). 

Technically speaking, one arrives at Eq. (jAip by using the background relations (fT^ ). (fTOh ) 
and (|21ll . together with the auxiliary linear evolution laws 

{daO^B) =-7HdaO^B ^ {da^'s) =-7HdaUJ% (^ 3 ) 
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and 


(d'^Aa)' =d^Aa-2Hd^Aa, 


(A4) 


which reflect the non-commutativity between convective and time derivatives (see expression m 
in footnote 3). Note that the last tw o relations stem f r om th e fact that {dbBa)' = —SHdbBa 
to first order (see also Appendix C in IDosodouIou et al (j2012l ll. Finally, we have obtained the 
homogeneous differential equation (jAip . by going back to Eq. (I27p and using it to eliminate the 
magnetic shear and vorticity terms. 


B. The case of a magnetic-free background 

So far, we have analysed the magnetic effects on the evolution of linear density perturbation, 
assuming a (Newtonian) Einstein-de Sitter background that contains a weak, random (i.e. statisti¬ 
cally homogeneous and isotropic) B-field. Here, primarily for comparison reasons, we will consider 
a pure (i.e. a magnetic-free) background cosmology. In the absence of the field, the zero-order equa¬ 
tions are given by the set (fT9h l- (fT9b l. Then, treating the magnetic pressure (H^) as a first-order 
quantity, the key linear variables are 

Aa = ^ daP , Za = adaO and = ” daB '^, (Bl) 

P P 

where Ba has now been normalised using the background density of the matter (compare to defi¬ 
nition p3b l in § 13.1|) . For simplicity, let us set the fluid pressure to zero and ignore the magnetic 
shear and vorticity terms. Under these assumptions, the above defined variables evolve according 
to the linear formulae 


Aa 

— —Za , 

Za = -2HZa-^KpAa-^d^Ba 

(B2) 

and 


Ba = -HBa , 

(B3) 

respectively. Proceeding in a way analogous to that outlined in Appendix A, the above expressions 
combine to the differential equation 


CO 

> 

T 

1 

10 d2A(„) 2 dA(„) 

(B4) 


dt^ 

3t dt^ dt ’ 

which accepts the solution 

A(n) ^ 

= Ci f^+Cs. 

(B5) 


The main feature here is the absence of explicit magnetic terms on the right-hand side (compare 
to solutions m, (El and ()A2I) of the weakly magnetised background). Instead, the H-field has 
simply added a (physically unimportant) constant mode in between the growing and the decaying 
modes of the non-magnetised case. Therefore, removing the magnetic field from the background 
model has weakened its overall effect, something clearly reflected in the above solution. Put another 
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way, incorporating a weak (fully random) i?-field in the background enhances its input by a small 
amount, which however is enough to allow for explicit magnetic terms in the solutions (see (j3ip 
and (|33|) b The latter clarify qualitatively the role of the field and also help to quantify its impact. 
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